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Motivated by the recent heat transport experiments in 2D antiferromagnets, such as La 2 CuC> 4 , 
where the exchange coupling J is larger than the Debye energy ©d, we discuss different types of 
relaxation processes for magnon heat current with a particular focus on coupling to 3D phonons. 

We study thermal conductivity by these in-plane magnetic excitations using two distinct techniques, 
Boltzmann formalism within the relaxation-time approximation and memory-function approach. 

Within these approaches, a close consideration is given to the scattering of magnons by both acoustic 
and optical branches of phonons. A remarkable accord between the two methods with regards to 
the asymptotic behavior of the effective relaxation rates is demonstrated. Additional scattering 
mechanisms, due to grain boundaries, impurities, and finite correlation length in the paramagnetic 
phase, are discussed and included in the calculations of the thermal conductivity k(T). Again, 
we demonstrate a close similarity of the results from the two techniques of calculating k(T). Our 
complementary approach strongly suggests that scattering from optical or zone-boundary phonons 
is important for magnon heat current relaxation in a high temperature window of Od < T J. 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.50.Ee, 72.20.Pa, 75.40.Gb 


I. INTRODUCTION 

After almost three decades of intensive studies, 
cuprates continue to attract significant interest because 
of their outstanding properties and due to the continued 
research effort in high-temperature superconductivity. 1 ' 
In particular, their magnetic properties remain inspira¬ 
tional, both as potentially responsible for the mechanism 
of superconductivitj® and on their own right as relevant 
to a larger class of low-dimensional antiferro ma gnets and 
as a test case for various theoretical models.- ^ However, 
understanding of some of the properties of the magnetic 
excitations in layered cuprates remains incomplete. This 
concerns interactions of such excitations with themselves 
and various other perturbations, as well as the role of 
such interactions in observable spectroscopic and trans¬ 
port phenomena. 

Earlier studies of thermal conductivity in La 2 CuOj^^ 
have demonstrated large contribution of magnetic exci¬ 
tations to the in-plane thermal transport, thus offering a 
unique window into their properties which are not easily 
accessible by other methods. More recent experimental 
advanced 1 -®! call for a deeper theoretical insight into the 
mechanisms of magnon heat current dissipation. This in¬ 
terest goes beyond a particular material and highlights 
a broader importance of general understanding of the 
transport phenomena in a wider class of antiferromag¬ 
nets. While scattering of magnons among themselves and 
due to fluctuations of the order-parameter in the para¬ 
magnetic state has received significant attention in the 
pastp^l the impact of phonons on magnon lifetime and 
other properties of magnetic excitations has only recently 
begun receiving attention: 13 

In this context, the current work focuses on the physics 
of magnon scattering and its role in magnetic thermal 


transport of quantum antiferromagnets. While our ap¬ 
proaches are generic, our analysis is strongly motivated 
by La 2 CuC >4 and related cuprates. Therefore, we con¬ 
centrate on the case of 2D magnons and superexchange 
coupling J large compared to the Debye energy Qjj, for 
temperatures T < 0.4J, where the validity of a magnon 
description is well-establishecP^. Moreover, most of our 
work is devoted to the magnon-phonon scattering to clar¬ 
ify the role of this less-studied relaxation mechanism. 
Additional effects, such as grain boundary and impurity 
scattering, and the role of finite correlation length are 
discussed less extensively. 

In addition to investigation of the transport proper¬ 
ties of layered cuprates, we also contribute to the formal 
development of transport theory by contrasting the re¬ 
sults from two complementary methods, the Boltzmann 
theory and the memory function technique. For the 
relaxation-time approximation within the Boltzmann ap¬ 
proach, we take advantage of the large energy scale of 
magnetic excitations and thus operate with asymptotic, 
long-wavelength expressions augmented by appropriate 
cut-offs and some well-justified modeling of the optical 
phonon spectra. For the magnon-phonon scattering we 
advocate the use of a simplified “effective phonon DoS” 
approach, which allows for straightforward yet fairly re¬ 
alistic calculations. In the memory-function approach, 
on the other hand, we maintain microscopic expressions 
for the magnon spectra valid in the entire Brillouin zone, 
while using coupling to a single dispersive phonon branch 
with a model dispersion, which introduces scattering 
on both acoustic and optical-like zone-boundary modes. 
The final results from both approaches, i.e. the thermal 
conductivities vs. temperature, are found to be remark¬ 
ably similar. Moreover, a complete agreement between 
both approaches on the power-law asymptotic regimes, 
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controlled by the optical and acoustic phonons, is also 
demonstrated. 

The paper is organized as follows: we begin with a gen¬ 
eral discussion of the qualitative features of the rnagnon- 
phonon scattering in Sec. |TT| followed by intuitively clear 
details of Boltzmann approach and calculations of k(T) 
within it in Sec. |III| We continue with the exposition 
of the memory-function approach and its results for the 
effective relaxation rates and thermal conductivity in 
Sec. [TV} A brief discussion of other scattering mecha¬ 
nisms and of their respective roles is given in the cor¬ 
responding sections devoted to the thermal conductivity 
calculations. Technical details, discussion of the physi¬ 
cal range of spin-phonon coupling, etc., are delegated to 
several Appendices. 


II. MODEL AND QUALITATIVE 
CONSIDERATIONS 

Generally speaking, a spin system on a lattice can al¬ 
ways be described by a Hamiltonian consisting of spin- 
only and lattice-only parts, TL S and TLph respectively , 
in addition to a coupling between them, which we will 
assume to be of magnetoelastic nature Hs-ph 

'H = 'H s -\- H p h + W-s- P h = r Ho + ?4-pii ■ (1) 

Having in mind La 2 Cu 04 and related cuprates, we take 
a simple, nearest-neighbor-dominated Heisenberg model 
with the superexchange constant J on a square lattice 
to be a faithful description of the 2D antiferromagnet of 
interest. Appendix[A]offers the standard linear spin-wave 
treatment of it leading to the free-magnon Hamiltonian 

'Us El > (2) 

k 

where £k is the magnon energy and h=kB = 1 from now 
on. 

The full phonon spectrum of La 2 Cu 04 is representa¬ 
tive of that of the other cuprates and comprises three 
acoustic and eighteen optical modes. With a typical 
bandwidth of each branch from 50K to 400K, together 
they span the range of energies reaching 900KP^ While 
in what follows we will model them in a simplified fash¬ 
ion, one can write their Hamiltonian as 

^P h ^ El CJ q/ a qf a qf > (3) 

where £ numerates branches of phonon excitations. 

A straightforward derivation yields the lowest-order 
magnon-phonon coupling in the following general form 
(see Appendix [B| for details) 

^s-ph = E E {t'k.k'.q&Ac (4) 

k.k'.q i 

+ 2 KEq (&PL + H.C.) } + a_ q <) , 
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FIG. 1: Magnon-phonon scattering diagrams. Solid lines are 
magnons, wavy lines are phonons, (a) Decay (phonon emis¬ 
sion), (b) recombination (phonon absorption), (c) anomalous 
process (absorption of magnon). 


where V kk , q and U kk '/ q are the “normal” and “anoma¬ 
lous” spin-phonon coupling vertices. For the coupling to 
acoustic and optical (or zone-boundary) phonons, they 
assume different asymptotic forms, discussed in the next 
Section and in Appendix |B| 

Since we are interested in the thermal transport by 
magnons, the following generic consideration is useful. 
Because the phonon Debye energy (0 D ss 400K) is much 
smaller than the magnon bandwidth (2.2J > 3000K), 
phonons can be assumed to be in thermal equilibrium, 
i.e., playing the role of a “bath”. This is well-justified for 
temperatures comparable to or above half of the phonon 
Debye energy T > 0d/ 2, which is about 200K in most 
cuprates. Then neither the momentum nor the energy 
of a magnon is conserved in the processes of magnon- 
phonon scattering, or, in other words, the momentum 
and energy are transferred from the magnon flow to the 
phonon bath. In that case, magnon relaxation time and 
transport times can be treated as the same. 

Thus, in contrast to the inter-magnon scattering, the 
magnon-phonon scattering channel is free from many re¬ 
strictions of the former and does not have the limitations 
on the phase space inherent to an Umklapp scenario. 
While the spin-lattice coupling constant may be small, 
phonons at temperatures T > 200K are abundant in the 
materials of interest. 

Another qualitative consideration, elaborated on in 
Sec. |III[ is that magnetic excitations are confined to lower 
dimensions than phonons (i.e., 2D vs 3D). In this case, 
momentum of phonons perpendicular to the 2D planes is 
not conserved, which also leads to fewer restrictions on 
the kinematics of the magnon-phonon scattering. 


III. BOLTZMANN APPROACH 

A. Relaxation rates 


Within the Boltzmann approach, the key element is the 
calculation of the scattering rates. As is argued above, 
transport and quasiparticle relaxation rates of magnons 
due to scattering on phonons should be the same. Then, 
in the lowest order of spin-phonon coupling, the problem 
is reduced to evaluation of the diagrams in Fig. [TJ which 
yield magnon relaxation rate by the standard diagram- 
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matic method 

^ = *£££ 


Tk 


(5) 


k' qy qx 


^ H k.k'.qj” (^q T ^-k / T 1) ^'k.k'+qy & (^k &k' ^q) 

T | hk' ,k,q | (n q ^k') 3 k',k+qy & (^k* ^k ^q) 

T | ^ k,k',q I (^k' ^q ) ^k'+k.qy & (-k' T £k ^q) ^ i 



with the first term corresponding to the diagram in 
Fig- 0a), in which phonon is emitted, and the second 
to Fig. 0b), in which phonon is absorbed. The third 
term is due to anomalous process, Fig. 0c), in which two 
magnons are absorbed and the phonon is emitted. We 
drop summation over the phonon branch index £ here, 
thus considering one branch of phonons at a time. 

The magnon-phonon vertices in 0 are the same as 
in Eq. (0 and for the first two terms they are related 
by a permutation of the initial and final states of the 
magnon. Note that the 2D momentum conservation in 
Eq- §, explicated by the <5 i,2+3’s, concerns only the 
in-plane momentum of the phonon q||, while the compo¬ 
nent perpendicular to the plane, qu, is not conserved. 
This is natural as magnons have infinite mass in the q^ 
direction. This feature is important for the future con¬ 
sideration and we separate sums over the components of 
phonon momenta q = (qy, gj_) in Eq. 0 explicitly. 


B. Approximations 

There are two approximations for La 2 CuC >4 and re¬ 
lated cuprates that follow from the fact that all relevant 
energies, T and ©d, are much smaller than the magnon 
bandwidth, W ~ 2.2 J: (i) magnon energies can be lin¬ 
earized 


£ k ~ v|k|, 


( 6 ) 


because for all practical purposes T <C J and hence |k| ~ 
T/J <C 1, (ii) similarly, magnon-phonon vertices for the 
optical and acoustic phonons are (F£ d k ~ Vk,k'. q ) 


Vk.k'.q « 

Fk.k'.q « • -J^=, 

vlql 


(7) 

( 8 ) 


where v is the magnon velocity (v = 1.15S\/2 J, lattice 
constant a = 1), and both g°. g* and are 0( J) < J (see 
Sec. IIIF and Appendices |B|and [c| for details)For the 
acoustic case, the phonon dispersion in the vertex is also 
linearized, in line with the Debye approximation. 

While magnon-phonon vertices in 0 and 0 can be 
proposed on general grounds, they can also be derived 
from realistic microscopic models of spin-phonon cou¬ 
pling, an exercise deferred to Appendix 0] The same 


FIG. 2: (Color online) Qualitative sketch of the magnon 
(£ k = ti|kj) and acoustic and optical phonon dispersions (o> q ’s) 
in the limit v^>c. Magnon and acoustic phonon energies are 
linearized. Schematics of the magnon [momentum k] absorb¬ 
ing a phonon [in-plane momentum qy] is also shown. 

Appendix also deals with the role of polarization of the 
3D phonons in the coupling to spins. The asymptotic 
k, k',q —> 0 form of these microscopic vertices agrees 
with 0 and 0, aside from some additional angular de¬ 
pendence that does not affect the results. The typical 
magnitude of coupling constants g sp ’s will be discussed 
in Sec. IIII FI 

We also make a note that coupling to the optical 
mode should typically involve a large in-plane momen¬ 
tum transfer qy, in which case a zone-boundary phonon 
from the nominally acoustic phonon branch is equivalent 
to the optical phonon. Such processes result in the scat¬ 
tering of a magnon from, e.g., the branch near k —> 0 to 
another branch near the AF ordering vector k' —> Qaf 
[= (7r,7r)]. Thus, in the following we do not distinguish 
between the optical and the zone-boundary phonons. 

The next line of approximations requires some qual¬ 
itative kinematic consideration. The linearization of 
magnon energies in 0-0 is well-justified for a typical 
£k ~T and already allows for some simplification in 0. 
Naturally, the “typical” range of momenta of magnons 
involved in the heat transport at relevant temperatures 
is limited by cxlT/v, which is <C 7r and concerns a small 
fraction of the Brillouin zone. Then, the typical in-plane 
component of the phonon momenta ]qii | = | ± k k'| for 
the phonons involved in the magnon scatterings in Fig. [l] 
must also be limited to the same range, |qy| <T/v, see 
the sketch in Fig. 0] Note, that if (acoustic) phonons 
would also be confined to 2D, this would imply that their 
typical energy needs to be much less than the energy of 
magnons: w q «c|q|| | ~ (c/»)TcT. 

However, for the 3D phonons the situation is radi¬ 
cally different. The energy conservations in Eq. 0 im¬ 
ply u q = |u(|k| ± |k'|)|, so for the typical magnon en¬ 
ergy T, the typical energy of a phonon is also T. This, 
in turn, implies that the out-of-plane component of the 
phonon momentum is much larger than the in-plane one, 
q_L>|q|||. This is particularly easy to see for the acous¬ 
tic phonon, for which the combination of the (in-plane) 
momentum and energy conservations yieldif 1 ^ 

c 2 (|k ± k'l 2 + ql) = V 2 (|k| ± |k '|) 2 , 


(9) 
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so for |k|, |k'|, |q|| | ~ T/v it follows that q± ~ T/c 
|q|||. For the optical phonon, w q « wo + a(q| + q±), the 
argument is simply that the typical |qn|~r/t;<Cl while 
q l does not have such restrictions. Note that for the 
optical phonon qy = qy ±Qaf is shifted by the in-plane 
AF-ordering vector. 

Altogether, this begs for the following approximation 


Wn 




9-L> 


( 10 ) 


simply neglecting the dependence of the phonon energy 
on the (small) in-plane momentum transfer |q|||. This 
approximation immediately simplifies Eq. © as the in¬ 
tegral over qy simply removes the in-plane momentum 
delta-functions, while the rest of the expression is inde¬ 
pendent of it (see a slightly more involved treatment of 
the case of the acoustic phonon later). Then the inte¬ 
gration over k' in ([ 5 ]) is simply removed by the energy 
conservation using linearized magnon energies. Lastly, 
the remaining integration over q± can be rewritten by 
introducing an effective DoS for phonons. Thus, one of 
the main results of this work is the development of an 
“effective” approach in which k- and T-dependencies of 
the relaxation rates in Eq. ([5]) are given by simple ID in¬ 
tegrals. A very close precision of this approach is demon¬ 
strated in Appendix [D] by a comparison with the direct 
integration in Eq. (|5j) without approximation of Eq. (10). 


C. Effective phonon DoS 


After using Eq. (10), the dependence of the integrand 
in Eq. ([5]) on q± is only through u q L , so it is natural to 
introduce an “effective” density-of-states of phonons 


D± (u) = J2 s ^~ u} ^'>- 


( 11 ) 


q± 


We would like to clarify that the “effective” DoS is not 
the full phonon DoS, but a ID version of it, which corre¬ 
sponds to the DoS of phonons with the vanishing in-plane 
momentum. 


1. Optical 


Given that the phonon spectrum of cuprates has more 
than a dosen of optical modes, covering the range from 
100K to 900Kp- 5 we reserve the right to model them in a 
more straightforward fashion. A sketch of such “model” 
densities of states is shown in Fig. [3] First model, which 
will be referred to as “Model I”, is just a constant DoS 
with the gap that corresponds to the lowest optical mode 


Dt M 


Q(uj - uj q) 

^max ^0 


( 12 ) 


The “Model II”, Fig. 131(b). corresponds to the optical 
mode in the sketch in Fig.Dlwith dispersion w q =wo+aq 2 






1 

1 

1 

1 

0)0 

0)0 

0)0 


(a) (b) (c) 


FIG. 3: (Color online) Qualitative pictu res of the effective 
DoS for the optical phonons, D ± (uj), Eq. (111. (a) “Model I” 
is the constant with the minimal energy cuo, (b) “Model II”, 
same with a more realistic square-root singularity at cuo, (c) 
“Model III” is for the “flat” mode at uio, Eqs. ( |12| ), (131, and 
(14), respectively. 


and includes a more realistic square-root singularity at ojq 
due to the ID nature of the “effective DoS” in (11) 


AiH = 


Q (cu - cu 0 ) 

27 xyj a{u — wq ) 


(13) 


with a= (w max — ujo)/tt 2 from normalization. “Model III” 


£>iiM =<5(w-w 0 ), 


(14) 


corresponds to the “flat” optical mode. We would like 
to note that the results for the thermal conductivity dis¬ 
cussed later are remarkably insensitive to the choice of 
the specific model for the phonon DoS. 


2. Acoustic 


The effective phonon DoS for the dispersion ui q± = cq± 


is a constant, same as the “Model I” in (12), but with no 
gap and an upper cutoff being 0£> (~ c) 


1 


£>acM = — 1 0(0D -w). 


(15) 


The validity of (151 is also restricted from below by 
w m i n = (c/v)e k at which q± sa |q|||, as we will discuss 
later. 


D. Scattering on the optical phonon 


The case of optical phonon is straightforward since 
magnon-phonon coupling within the approximation 0 
does not depend on the phonon momentum. Using lin¬ 
earized form of magnon dispersion and magnon-phonon 


vertex from © and f7]), approximation of Eq. ( fIo| ), and 
replacing (>0 k , with k!dk! which takes into account 


two magnon modes per Brillouin zone of the square lat- 
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tice, yields 


1 

7k 


<“> 

X |©(ek - w) (£k - w) 2 (n(w) + n(e k - w) + 1^ 

+ (e k + w) 2 (n(u) - n(e k + w)) 

+ 0(w - e k ) (w - £k ) 2 (n{u - £k) - n(w)) |, 


where the first term is from the phonon-emission dia¬ 
gram in Fig.[l|a), the second is due to phonon-absorption 
in Fig. 0b) , and the third is the anomalous term in 
Fig. 0 c), first, second, and third terms in Eq. ([5|, respec¬ 
tively. For the phonon-emission term the phonon energy 
is limited from above by the energy of the magnon that 
emits it, £k, and in the anomalous term the situation is 
reversed as the phonon energy must exceed that of the 
magnon, so the integration is limited from below. In 
(16), uj is also restricted implicitly through the DoS by 


wo, the lowest energy of the optical mode, and by w max , 
the highest energy of phonon bands. It is assumed that 
£k < w max - Thus, after all the legitimate approximations 
discussed above are implemented, we have a compact ex¬ 
pression of the relaxation rate in terms of ID frequency 
integrals, Eq. (16). Note that for the “Model III” (flat 


phonon mode) of the effective phonon DoS this integral 
in Eq. © is trivially removed and the relaxation rate 
is given by a compact analytical expression, presented in 
Appendix [D] 

In the limit of low temperatures, T, £k<SCwo, phonon- 


emission term in 1/r in (161 is strictly zero as the magnon 


with £k < wo cannot emit an optical phonon, and the two 
remaining terms are exponentially small, 

For higher temperatures T > ©£>, using the hierarchy 
of scales T > £k > wo (see Ref. USD and n(w)~T/w yields 
the asymptotics for the first two terms in 


Ai) 


oc 


Tel 


OC 


T 2 4 


, k V 2 UJ 0 ’ T] ( 2 ) ~ U 2 W 0 


(17) 


where we used g°^ ~ v. This is valid for Model I and 
Model II of the effective phonon DoS discussed above, 
while for the Model III (flat phonon), l/M 2 ) has the same 
asymptotics l/r^b Although it is natural to expect that 
the heat is conducted largely by thermalized magnons 
with a “typical” £k ~T, this is not exactly so in our case, 
because the distribution of the heat conducting magnons 
“leans” towards lower energies as we shall discuss later. 
Nevertheless, it is worth pointing out that “on a ther¬ 
mal shell”, i.e., at £k = T, where the memory-function 
and Boltzmann approaches can be consistently compared 
with each other, both expressions in © yield the same 


1 T 4 

-r CX 77- . 

r °pt V 2 U 1 0 


(18) 



FIG. 4: (Color online) T-dependence o f th e magnon relax¬ 
ation rate on the optical phonons, Eq. (161, for e^ = T and 


using effective phonon DoS models in Eqs. (12 I (141. The 


results are nor ma lized to the high-temperature asymptotic 
behavior, Eq. (181, [T ,4 /v 2 wo]. Parameters are as discussed 
in text and indicated in the graph, w max = 950K for Models I 
and II. The vertical axis is in units of (g° p 1 /v) 2 . Inset shows 


individual contributions of the three terms in ( |D1| ) [diagrams 
in Fig. [I] for Model III. The results for the relaxation rate 
due to acoustic phonons, Eq. (211, are also shown for a repre¬ 


sentative choice of the Debye energy 0£>=4OOK. J=1500K. 


This coincides with the results of the memory-function 
approach discussed later. We note that the phonon- 
emission is always subleading to the absorption term, 18 
and the anomalous term is negligible at high T. 

Our Figs. [4] and [14] demonstrate these asymptotic 
trends explicitly. They also show a close quantitative 
similarity of the magnon relaxation rates obtained from 
Eq. (16) using three different models for the effective 
phonon DoS, Eqs. (12)-(14). In Fig. [4j the T-dependence 
is shown for the relaxation rate for £k = T, i.e., on the 
thermal shell. The results are normalized to T 4 /'c 2 wo to 
make the high-tenrperature asymptotic behavior of (18) 


apparent. The vertical axis is in units of (g°^/v) 2 . In¬ 
set shows individual contributions of the three terms in 


(16) for the “Model III” (14). For this model, the con¬ 


tributions of the phonon-emission and anomalous terms 
[first and third in (D1)] are explicitly limited by the step- 
functions, resulting in a kink at T = wo, while for the 
other models these are smoothed out. For the sake of a 
comparison with the acoustic phonons in the next Sec¬ 
tion, we also note that the results for the scattering on 
the optical phonon in Figs. |4| and [14| likely underestimate 
the effect by a numerical factor about 2, due to scattering 
of magnons between k —»• 0 and k —> Q af branches. 


E. Acoustic phonon 

Below we demonstrate that the effective phonon DoS 
approach can be successfully extended to the considera¬ 
tion of magnon scattering off the acoustic phonons. 
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The crucial difference of the acoustic phonon scatter¬ 
ing is the form of the magnon-phonon coupling in Eq. ([8]) , 
which leads to an extra factor |q| | 2 /|q| in the scattering 
probability compared with the optical phonon case. A 
naive power-counting, together with the kinematic con- 
IIIB suggest that the typical |qy | ~T/u 
(u/c)|qii|, translating this extra factor 


sideration in Sec. 
while |q| « <z_l « 


into cTjv for the relaxation rate. Using that the phonon 
sound velocity c~0b~uo, this would yield the “thermal 
shell” estimate for the relaxation rate 1 /T ac ~ T 5 /v 4 , to 
be contrasted with the asymptotic expression for the op¬ 
tical case (18). The situation is more delicate, however, 


as the scattering in the present case is, in fact, dominated 
by the low-w phonons. Technically, the integral over the 
phonon energies diverges as 1/w and must be cut off at 
Wmin = (c/v)e k, where q± becomes ~ |q|||. Altogether, 
assuming g *this gives the “thermal shell” estimate 
for the acoustic case as 


1 


oc 


Qdv 3 ’ 


(19) 


which is valid for both T> 0£> and T < ©£>. Note, that 


the result (19) is the same as in the memory-function 


consideration, Sec. [IV 

A rigorous derivation of this result from Eq. ([5]) needs 
a slightly more delicate treatment of the magnon-phonon 
coupling. We use |qn | 2 = |k' —k| 2 and, according to ap¬ 
proximation in Eq. (10), |q|~< 7 _i_. Then the extra factor 


|q| in the scattering probability reads 


,'|2 


|k| 2 — 2|k'||k| cost/? 
Q± 


(?) 


st + 


u q± 


, ( 20 ) 


where ip is the angle between k' and k and, because of the 
approximation of Eq. (10), the term with cos</3 averages 


to zero upon the integration over this angle. With the 
result in Eq. (201 and using the effective DoS model for 


the acoustic branch introduced in Eq. (15), we rewrite 
the relaxation rate in ([5| as 

x ((e k - w ) 4 + el (e k - w) 2 ) (n(u) + n(e k - w) + l) 

+ ((e k + w ) 4 + el (e k + w) 2 ) (n(u) - n(e k + w)) 

+0(o; — e k ) 

x ((w - e k ) 4 + el (w - e k ) 2 ^ (n(u - e k ) - n(w)) |, 

where the three terms are the phonon-emission, phonon- 
absorption, and anomalous terms in Fig. [l|a)-(c) and in 
Eq. (|5|, and w m i n = (c/u)£ k as before. 

We hrst note that the contribution of the anomalous 


term in the acoustic phonon case of Eq. (21) is by a factor 


of c/v smaller than that of the other two terms. The sub¬ 
tle reason for that is in the threshold nature of the pro¬ 
cess: the lowest possible phonon energy is e k , not w m i n , 


which gives the thermal-shell estimate I/t® ~ T 5 /u 4 , 


much less than the result in (19). We, therefore, give the 


asymptotic consideration only to the first two terms. 


Because the integral over the phonon energy in (21) 


is infrared-divergent and thus is dominated by the low- 
energy phonons, the asymptotic behavior of both terms 
is the same at low and high temperatures (T < 0£> and 
T > 0_d). For e k being of the same order as T, both 
terms in (21) yield an estimate of the relaxation rate 


O). (22) 


1 

d2) 


5sp 

V 


Tei 


which is in accord with the thermal-shell answer (19). 


As we discuss below, the phonon-absorption term has a 
more complicated e k -dependence for e k <CT. 


The asymptotic result (19) should be compared with 
the high-temperature (T > ojq ~ 0d) estimate for the 


optical phonon case, Eq. (18). The ratio of (19) to the 


latter is T/v, which should imply that the contribution of 
the scattering on acoustic phonons is a relatively minor 
effect in this temperature regime. A direct comparison is 
provided in Fig. |4j where the black line is obtained from 
(21), without the use of the asymptotics. This line clearly 


indicates that the asymptotic consideration of Eqs. (19) 


and (22) is correct and that the relaxation rate on ther- 
malized acoustic phonons follows T 5 power law. 

On a closer inspection of Fig. [4] we should note, first, 
that the dominance of the acoustic phonon scattering 
in the low-T regime concerns a really small region of 
T<0£)/4 and is unlikely to be seen in the thermal con¬ 
ductivity of La 2 Cu 04 as this regime is known to be dom¬ 
inated by the grain-boundary scattering. 9 Second, at the 
higher T, the relaxation rate by acoustic phonons seems 
to exceed the one by optical phonons, at least for some 
of the models of their DoS. This is likely to be due to a 
neglect of the numerical factor difference in the coupling 
strength to acoustic and optical phonons in our effective 
vertices Q and 0, and an underestimate of the optical 
phonon scattering rate by a numerical factor mentioned 
in Sec. IIIIDI 


F. Summary of the phonon-scattering mechanisms 

Here we would like to summarize our considerations 
of the magnon-phonon scattering and to take a broader 
view of its implication for the thermal conductivity. 


1. Smallness ofg’s 

First, while the perturbative character of our treat¬ 
ment of magnon-phonon scattering is implied by the use 
of the lowest Born approximation in Fig. [l] and Eq. ([5|, 
we would like to make it explicit that the physical range 
of the phenomenological magnon-phonon constants we 
introduce in |7| and Q is g sp /v <C 1. At first glance 
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this may be surprising as the dependence of the su¬ 
perexchange constants on the interatomic distance is of¬ 
ten rather sharp and is governed by some high power 
of the distance, leading to estimates dJ/da ~ 7 J with 
7 ~ 10 — 20, see Ref. [TB] However, this largeness is offset 
by the smallness of a characteristic atomic displacement 
associated with phonons,^ Ify/mOp ~ 1/100, see Ap¬ 
pendices [B] and [C] on how the two factors appear together 
within a microscopic approach. 


modes at the temperatures T around wo,i=wo,i/2, as in 

Fig-El 

We remark that thermal conductivity obtained with 
this effective 1 /T k and appropriate set of parameters can 
be made virtually indistinguishable from the ones using 
more elaborate expressions from Secs. mu and |HIEl 

G. Other scattering mechanisms 


2. k-dependence 

Second, we note the importance of the ^"dependence 
of the relaxation time. For bosons with e k ss u|k|, one 
can estimate thermal conductivity 


r T/v 


K OC 


Tk dk. 


(23) 


According to the preceding sections, in magnon-phonon 
l/r k the lowest power is k 2 . In a similar situation in 
1D,E3 this leads to a strong infrared divergence of the 
spin component of the thermal conductivity. This means 
that the spin-phonon scattering is not sufficient to ren¬ 
der conductivity finite and one needs to take into con¬ 
sideration other scattering mechanisms. In 2D, the inte¬ 
gral in (23) still has a weak (logarithmic) divergence for 
r~/c -2 , so the (grain-)boundary scattering is sufficient 
to mitigate it. One of the implications of this is that the 
distribution of magnons that carry the heat most effec¬ 
tively is not centered at energies of order T, but is shifted 
toward lower energies. More importantly, this considera¬ 
tion means that with the e k -dependence of the relaxation 
rates obtained in Secs. |III D and III E[ magnon-phonon 
scattering cannot be the only scattering mechanism and 
thus must be accompanied by a boundary-like scattering 
in order to render magnon heat conductivity finite. 


3. Effective 1/r 


Lastly, given how close the results for different models 
of the phonon spectra conform to the asymptotic expres¬ 
sions for 1 /r k in 0 and | 22 | ) , it is tempting to introduce 
a simplified, “effective” expression for the magnon relax¬ 
ation rate on optical and acoustic phonons that contains 
a minimal number of parameters 


f \ ' opt 

/iff ~Z^ a i 


r 2 k 2 

0 (T — ojq 7 )- 1 - 

Wo ,j 


vTk 4 


(24) 


where a ac = {g^/v) 2 and a° pt = (g° p \/v) 2 are the dimen¬ 
sionless coupling constants to the acoustic and the zth 
optical mode which has the (lowest) energy wo,i- The 
“pseudo”-step-function ©(a;) is introduced to mimic— 
the exponential “turn-on” of the scattering on optical 


1. Grain-boundary 


As is discussed in previous section, boundary scatter¬ 
ing is essential to mitigate residual infrared divergence of 
magnon n m if only scattering on phonons is considered. 
One can expect that in a 2D magnetic lattice even rel¬ 
atively weak dislocation-like defects are likely to act as 
strong boundaries for magnon propagation, similarly to 
the effect of crystal grain boundaries on phonons. Experi¬ 
mentally, the grain-boundary scattering is known to dom¬ 
inate entirely the low-temperature (T < 200K) magnon 
thermal conductivity of La 2 Cu 04 P The corresponding 
relaxation rate is simply 


v 

V 


(25) 


where L is the characteristic size of the grain. The typ¬ 
ical grain size quoted in Ref. [5] is L ~ 150 lattice spac- 
ings and is in agreement with the other measurements in 
La 2 Cu 04 . 


2. Correlation length 


In the paramagnetic state above the Neel temperature, 
which is non-zero in most unfrustrated 2D AFs because 
of the interplane interactions and/or small anisotropies, 
finite spin-spin correlation length can be expected to rep¬ 
resent a natural “cut-off” boundary for magnon propa¬ 
gation because magnons are the spin-flips in an ordered 
structure. This expectation is in agreement with a num¬ 
ber of studies which have pointed out that the dynam¬ 
ics of the antiferromagnets is fully diffusive at scales be¬ 
yond the correlation length £(T), i.e., that the propa¬ 
gating magnons are overdamped by fluctuations at dis¬ 
tances A>£(T).™ai n effect, the correlation length acts 
a temperature-dependent size of an order-parameter do¬ 
main for magnon propagation. 

The exponential dependence of the 2D correlation 
length on J/T has been verified experimentally in 
La 2 Cu 04 and in the other 2 D antiferromagnetsP^H 25 ^ 
and is supported by extensive theoretical and numer¬ 
ical Quantum Monte-Carlo (QMC) calculations ! 21 ! 26 ! 27 ! 
An approximate analytical expression of it for the spin-1, 
nearest-neighbor Heisenberg antiferromagnet on a square 
lattice iJ3 




1.13J 
2.26 J + T 


e 113J/T , 


(26) 
























in units of lattice spacings. While it was argued 
that the characteristic length-scale that defines magnon 
lifetime should also contain a temperature-dependent 
prefactoi ^ 21 l 23 l 24 l 2 ^ as well as correctional factors,^ we 
simply suggest the following scattering rate by interpret¬ 
ing correlation length as a mean-free path 


1 v 

t^wy 


(27) 


These seemingly naive expression and expectation are 
supported by the studies of the T-dependence of 
the magnon linewidth in copper-formate-tetradeuterate 
(CFTD), another model spin-f square-lattice antiferro- 
magnet, by inelastic neutron scattering and QMCP^ In 
fact, the relaxation rate in the form of Eq. (27) has been 
suggested in Ref. [25} which has demonstrated that in 
both experiment and QMC results the magnon linewidth 
is well described by (27) at T>0.2J. 


Another effect of the finite correlation length is an ef¬ 
fective gap A;r = u/2£ in the magnon spectrum 14 


£k 


W* 2 +£- 2 cn/4, 


(28) 


which, obviously, affects the contribution of the long- 
wavelength excitations to thermal conductivity. 


3. Lattice disorder 


In a recent studyp^ the following scenario has been put 
forward. Since the ionic motion in the cuprates is much 
slower than the superexchange processes among spins, 
the former must cause a significant variation of couplings 
J due to zero-point and thermal lattice fluctuations. The 
effect was estimated from the high-resolution neutron 
diffraction and the zero-point motion was found respon¬ 
sible for the distribution of the width SJ ~0.1j] 13 In a 
sense, this implies that magnons propagate in a medium 
with the velocity that is randomly varying around the 
mean value with a distribution given by SJ(T). One can 
approximate this effect by the T-dependent static ran¬ 
dom lattice disorder with the lowest-order magnon re¬ 
laxation rate due to this mechanism given by 


1 

-T-lat 

T k 


SJ 2 (T) 7 3 


(29) 


where the temperature-dependent disorder strength can 
be modeled as SJ(T) = SJ(0)y/l + 4 T/<3 D to in¬ 
terpolate between the amplitude of zero-point and 
temperature-induced lattice fluctuations within the De¬ 
bye approximation. OIHD 

Related to this mechanism are two other possi¬ 
ble sources of magnon scattering at high tempera¬ 
tures that are harder to estimate. First is specific 
to La 2 Cu 04 , which exhibits orthorhombic-to-tetragonal 
structural phase transition at about 525K associated with 
a softening of a phonon mode? 1 This transition may have 


a direct impact on the values of superexchange constants 
and also enhance magnon-phonon scattering involving 
the mode that is being softened. However, it is hard 
to quantify both without a microscopic insight. 

Second is a significant decrease and eventual collapse 
of J at high enough temperatures, advocated in Ref. \W, 
as an ultimate result of the thermal expansion. With the 
large value of J and an apparent insignificant impact of 
the expansion on the average value of J in the 0 — 300K 
ranged in La 2 Cu 04 , it is hard to estimate at what T 
such dramatic effects can be expected to onset. 


4- Magnon-magnon scattering 


Last but not the least is the effect of magnon-magnon 
scattering. Since these require an explicit momentum 
dissipation to contribute to conductivities, the standard 
Umklapp process leading to a scattering of a typical low- 
energy magnon with the momentum k must involve a 
high-energy magnon with the energy £ max ~2.2 J. Then, 
one can suggest an ansatz 


1 

mm,U 

r k 


Jk e _2 2J / T . 


(30) 


This can be seen as an upper limit estimate for the 
standard Umklapp scattering rate as it neglects possi¬ 
ble smallness of the matrix element and only takes into 
account the smallness of k from the initial state. 

There exists a possibility of an unconventional “low- 
energy Umklapp”, because magnons in La 2 CuC >4 have 
two branches, at k = (0,0) and k = (7r, 7 t), so that some 
of the scattering between these branches may carry away 
large momenta. This logic seems to be implied in a re¬ 
cent calculation in Ref. [52] In that case, the correspond¬ 
ing relaxation rate can be expected to display a power- 
law behavior, similar to the “normal” magnon-magnon 
scattering 12 


T 2 k 


(31) 


which, on “thermal shell”, differs by a factor T/oj o from 
the magnon-phonon estimate in (18). This implies that 


for sufficiently high temperatures, T > magnon- 
phonon scattering must be more important than the 
inter-magnon scattering. 

Moreover, this expression must carry a small prefactor 
as it neglects the smallnesses of the fraction of the Umk¬ 
lapp vs normal processes and of the phase space for scat¬ 
tering. The apparent inability of the magnon-magnon 
scattering theory of Ref. [22] to fit the data for La 2 CuC >4 
from Ref. 2] beyond the boundary-controlled regime can 
be seen as an indirect evidence of the relative unimpor¬ 
tance of this type of scattering for the magnon thermal 
conductivity in large-J materials. We thus exclude it 
from the subsequent consideration. 
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FIG. 5; (Color on line ) Magnon relaxation rates from 

Eqs. |24| ), (251, (271, (29 1 , and (30) normalized to T 4 /v 2 Qd, 
Eq. (181, at £k = 4 T. Parameters are as described in text. 


5. Comparison 


To gain a qualitative expectation for the relative con¬ 
tribution of the discussed scattering mechanisms in dif¬ 
ferent temperature regimes, we compile the results from 
Eqs. (24), (25), (27), (29), and (30) in our Fig. [5j 
which shows relaxation rates normalized to the high- 
T asymptote of the optical-phonon scattering in (l8| ) 
[l/r opt =T 4 /'(; 2 0fl], as in Fig. 4j For the magnon-phonon 
scattering (24) we have dimensionless coupling constants 


opt 


= a ac = 0.1 and a 2 Pt = 0.3 with wo,i = B_d = 400K 
and wo ,2 = 900K. For the boundary scattering (25) the 
scattering length is L = 300 lattice spacings, and for the 
correlation length effect (27) we use the expression for 
£(T) inEq. (26). Temperature-dependent lattice disorder 
coupling in (29) is chosen with 5J{ 0) =0.1J to roughly 
match the results of Ref. in J = 1500K as before and 
magnon energy £k is chosen to be e k = 4T, the value 
of energy which roughly corresponds to the maximum of 
the thermal population of magnons at a given T. 


One can see in Fig. [5] that the scattering at low-T 
(T < 0 _d/ 2) is controlled entirely by the grain bound¬ 
aries. We note that the magnon-phonon scattering and 
the lattice-disorder scattering have a substantial depen¬ 
dence on £k and thus lead to a stronger scattering for 
magnons with £k > T, but to a weaker scattering for 
£k < T. Thus, in Fig. [5] at intermediate and high-T 
(T > 0 d / 2) the dominant contribution is due to finite 
correlation length and phonons. The magnon-phonon 
scattering is at least as strong as the scattering due to 
correlation length at intermediate T and even exceeds it 
for the higher T for the given choice of £k- While lattice 
disorder effect is also significant, it remains secondary 
and, given its nearly perfect T 4 asymptotic behavior, may 
be incorporated into the magnon-phonon coupling to an 
optical mode. Thus, while the lower-energy magnons are 
scattered almost exclusively by the boundaries and finite 
correlation length, at T>Qp> the higher-energy magnons 


are strongly scattered by phonons. 

Altogether, grain-boundary, correlation length, and 
magnon-phonon scatterings are the leading mechanisms 
of magnon energy relaxation and, therefore, have to be 
included in the consideration of thermal conductivity of¬ 
fered next. 


H. Thermal conductivity 


Within the Boltzmann formalism thermal conductivity 
by bosonic excitations in 2D is 


k 


(cos^Uk) 2 ( y ) («k + «k)T k , (32) 


where <p is the angle between Vk and the current. Using 
the linearized form of magnon dispersion © and the fact 
that all of the discussed relaxation rates are isotropic in 
k we can simplify (32) to a ID integral 


'2 _i 2 /*^max 

—— / dx 
27r J n 


z 3 e x 
(e x - iy 


t(x,T), 


(33) 


where x = s^/T and x max = £ max /T with £ max = vV2rr, in 
which we used the Debye-like approximation for magnons 
and also accounted for two magnon modes in the Bril- 
louin zone. 

At finite correlation leng th £, the magnon spectrum 
opens a gap^ according to (28), which results in a mod¬ 


ification of the expression for the thermal conductivity 

rj^2 /*3?max / ^,2 


2tt 


■i'H-lpf) (34) 


where x m - ln = A^/T with A^=v/2 £(T). 


1. Comparison 


Finally, we present the results of our calculations 
of magnon thermal conductivity in Figs. [6] and [TJ in 
which we demonstrate the effects of boundary scattering 
(25) together with the finite correlation length and/or 
magnon-phonon relaxation mechanisms. In Fig. [6j the 
upper solid curve shows k(T) for the case when the scat¬ 
tering is only by the grain boundaries (251 with L = 300 
and by phonons (|24|) with the same coupling parame¬ 


ters as in Fig. [5j The middle solid curve is for n(T) 
due to grain boundaries (25) and the correlation length 


(27). The corresponding dashed lines are for the same 


cases, but with an additional effect of the finite-T gap in 
the magnon spectrum due to finite correlation length, 
Eq. (341. The effect is minimal on the middle curve, but 


is rather dramatic in the phonon-scattering case. This 
is due to strong £k-dependence of the magnon-phonon 
scattering discussed in Sec. |III G~1 which changes substan¬ 
tially the magnon population contributing to the heat 
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FIG. 6: (Color online) Magnon k(T) for the scattering by 
the grain boundaries, L = 300, and either the correlation 
length (middle c urve s), Eq. (271, or phonon scattering (up¬ 
per curves), Eq. (24 1 . Lower curves combine the effects of all 


three, grain-boundary, phonon, and correlation length scat¬ 
tering mechanisms. Parameters for the phonon scattering are 
as in Fig. [5] D ashed lines show the effect of the gap in magnon 
spectrum (|34|). 


current. Namely, in the “boundary+£(T)” case, the typ¬ 
ical magnon in (341 is a thermalized one, £k ~ T, so 


the exponentially small gap Af(T) of (28) does not af¬ 
fect it. In the “boundary+phonon” caseTThe high-energy 
magnons are scattered strongly by phonons, see Fig. [5j 
while the heat-conducting population of magnons leans 
strongly to the low energies, hence a dramatic impact of 
opening the gap. 

Lower curves shows the result of combining all three, 
the grain-boundary, the phonon, and the correlation 
length scattering mechanisms. Clearly, the phonon re¬ 
laxation mechanism leads to a stronger scattering at the 
higher-£k part of the heat-conducting magnon popula¬ 
tion, reducing the overall conductivity, while the longer- 
wavelength part is now controlled by the correlation 
length, evidenced by a weak sensitivity of k(T) to the 
gap in the magnon spectrum (dashed curves). 

Fig. [^complements this study with the consideration of 
three different scenarios: the magnon-phonon coupling is 
to (i) two optical modes, one at wo,i =400K and the other 
at w 0 ,2 = 900K (an analog of the high-energy stretching 
mode), (ii) to an acoustic branch of phonons only, and 
(iii) to both (same curve as in Fig. [6]). One can see that 
the overall effect on k(T) is very similar in all three cases. 
This is natural in the light of the preceding discussion 
that points to the significant scattering effect by phonons 
at higher energies, which can be expected to be compara¬ 
ble regardless of the nature of the phonon branches. We 
would also like to point out that we have performed the 
same calculation of k(T) in (34) using not the “effective” 
magnon-phonon relaxation rate, proposed in Sec. IIIF 
in Eq. (24), but the relaxation rates in Secs. HID and 


HIE Eqs. (16) and (21), using different effective phonon 


DoS models, previously matched to a direct numerical 


FIG. 7: (Color online) Magnon k(T) with all three key scat¬ 
terings, same as in Fig. [6] For the magnon-phonon scattering, 
results for coupling to only optical, only acoustic, and to both 
are shown. 

integration in Eq. ([5]). Aside for the adjustments of the 
magnon-phonon coupling constants needed in the cases 
when phonon spectral weight is spread over a substan¬ 
tial energy range, these calculations have produced the 
results that are virtually identical to the ones in Figs. [6] 
and0 

Altogether, magnon thermal conductivity in La 2 CuC >4 
should be largely controlled by the “boundary-like” scat¬ 
terings, coming from either the real grain-boundaries or 
the 2D correlation length, with a substantial correction 
from the magnon-phonon scattering, affecting k(T) at 
intermediate and high temperatures. 


2. Remarks 


For the latter regime, there are two additional remarks 
that need to be made. First, according to Eq. (26), at 


800K («0.5J) the correlation length is of order of 3 lat¬ 
tice spacings. Is the magnon picture still applicable and 
is it realistic to expect further reduction of the mean-free 
path by a phonon scattering? We believe that although 
the magnon description at T>0.5J may be an extrapola¬ 
tion, such a correlation length corresponds to a “patch” 
of about 30 spins, which is known to be well-described 
by the spin-waves in finite-cluster studies. A somewhat 
more relevant quantity is the effective gap in the magnon 
spectrum at this temperature, « J/4, which is still 
considerably smaller than the magnon bandwidth (2 J). 
Given that the magnon-phonon scattering mostly affects 
magnons with £k > T, it seems perfectly legitimate that 
the phonons are the source of a further shortening of the 
mean-free path in this regime. 

Second, the use of the bosonic description of the ther¬ 
modynamics of spin excitations of the S = \ Heisenberg 
model on a square lattice is limited by the temperatures 
of order 0.6,7, at which the specific heat shows a broad 
maximum 30 while within the bosonic description specific 
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heat saturates at a somewhat higher T. For a compar¬ 
ison of our results with experimental data this implies 
that the high-temperature tail of our k(T) should only 
serve as an upper limit estimate of the reality. 


IV. MEMORY FUNCTION APPROACH 

We now turn to the memory function approach, which 
does not proceed via one particle excitations, but focuses 
on the dynamics of the current directly, as we will sketch 
next. We start from the magnon heat current in terms 
of Bogoliubov quasiparticles 


j = 51 £kVk ^A - ( 35 ) 

k 

where Vk = Vk£k is the velocity of the magnon. Fol¬ 
lowing the memory function method^, the dynamical 
thermal conductivity tensor at frequency z = to + *0 + is 


VTk^ v = 



i 

z — C 




M (z)x~ 1 ' 


fiiy 


(36) 


where V is the volume, and C = £ s + £ p h + £s-ph = £o + 
£ s _ p h is the Liouville operator CA = [H, A], comprising 
spin, phonon, and spin-phonon parts, see (I1J). (A\B) = 
Jq ( A + (X)B)d\ — (3(A + )(B) is the Mori’s scalar product 
with A(A) = Ae~ xu and f3 = 1/T is the inverse tem¬ 
perature. x is the isothermal heat current susceptibility, 
Xnv = (ijulii/), and M^(z) = {Cjp\{z - QC)~ l QCj v ) is 
the memory function matrix. Q is a projector perpen¬ 
dicular to the heat current in terms of Mori’s product 

Q = 1 - 

We will evaluate |36| for ui —» 0 to order 0(H^_ ph ). 
Within perturbation theory 31 'to the leading order in the 
spin-phonon coupling, the memory matrix is given by 


— V [n s - ph ,j„]) , (37) 

^ Ao / o 

where the subscript 0 refers to Mori’s product and ther¬ 
mal averaging with respect to the canonical ensemble of 
the system at zero coupling, "Hg-ph = 0. This subscript 
will be dropped hereafter. Since M pv ~ 0(H% h ) al¬ 
ready, the static current susceptibility is needed only for 
Hs-ph = 0. With [Ho Jp] = 0 and (jp) = 0 

Xn» = ^r( kk) = 5I £ k v k^( £ k) (l + n(e k )) , 

(38) 

where n(s k) = l/(e £k//T — 1) is the Bose-function. 


Aifiv{z) — ( [H s — ph)J/i] 


A. Relaxation time approximation 

The relaxation-time approximation corresponds to re¬ 
placing M(z)/x by — i/r with a phenomenological scat¬ 


tering time r. In that case (36) reads 


= ^2 51 £ k v k n ( £ k) (1 + n{e k )) 


(39) 


with Kp U = This is identical to the standard result 

from kinetic theory. For T <C J, Eq. (39) yields kocT 2 , 


as expected for bosons in 2D. While the case T J is 
unphysical, we note that then k oc const. 


We note that (39) is completely identical to (32) for 


a momentum-independent scattering rate. In principle, 
the memory function approach in Eq. (36) can be for¬ 


mulated with a particle-hole type of observable, jk q = 
£kVk/3k/3 k + q ; in order to explicitly analyze momentum- 
dependent current relaxations rates. We will not pursue 
this direction here. 


B. Spin-phonon relaxation rate 


Here we would like to restrict ourselves to the case 
of acoustic phonons. This reduces the problem to a 
monoatomic Bravais lattice with ionic masses m where 
magneto elastic couplings results from stretching of the 
nearest-neighbor exchange bonds along the in-plane bond 
directions, the case considered in detail in Appendix |B| 
Lattice deformations perpendicular to the bond direction 
lead to higher-order couplings. 

First, we use explicit form of "H s _ p h in Q to find the 
commutator [H s -ph,ju\, which has a meaning of a force 

[fts-ph, 3»] = 2^ E B l F fc, («E + %e) , (40) 

kq^ 


where k' = k + qy imposed by the conservation of the in¬ 
plane component of the momentum, £ numerates acoustic 
phonon branches, £?£ = (/3k,/3 k )> and the matrix 


= 

£ k.q - 




V 


k,k+qi| ,q I 


(41) 


where we have introduced the shorthand notation for the 
energy current of a single magnon mode ek = £kVk and 
used inversion symmetry e_k = £k, v -k = — Vk- The 
2x2 “vertex matrix” V£ k , is built from the diagonal el¬ 
ements responsi ble fo r the “normal” scattering processes, 


Vjf.k'.q given in (|B5|), and the off-diagonal, “anomalous” 


terms V^f given in (jB6|. 

As usuaPS (37) can be rewritten as zMp V (z) = 


XpA z ) - X% where XpA z ) is the retarded dynam¬ 
ical force-force susceptibility, resulting from analytic 
continuation of the imaginary-time Greens function 
X% v (t) = (T t ( [Hs-ph, jp] + (t)['H s _ p h, jp\)) , and x pv = 


H s -ph, jn\\[H s -ph, jp]) is the isothermal force-force sus¬ 
ceptibility. Due to tetragonal symmetry, all of these 
quantities are diagonal with respect to /i, v. 

For thermal transport we focus on the DC limit w—»0. 
Since Re[M A1!/ (iO + )] = 0 and Imy^ = 0, we have to cal¬ 
culate Im [Xui/( z = ®0 + )] on ly- This can be done using the 
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lu~k 

3kq;n 


Jkq;o 


FIG. 8: (Color online) Force-force susceptibility diagrams. 
Solid lines are magnons, wavy lines are phonons, (a) nor¬ 
mal contribution (magnon emission-absorption), (b) anoma¬ 
lous contribution (two magnon emission (absorption)). Filled 
circles in (a) are normal and in (b) anomalous force-vertices. 


two diagrams in Fig. [8j which comprise normal magnon 
emission-absorption and anomalous two-magnon emis¬ 
sion (absorption) processes, i.e. x£™( z ) and Xp“ 1,2 ( z )- 
After some algebra we get 

M^(iO + ) = -i n(u4) (42) 

k,q,^ 

x(2|<J 2 ^(^k) (l - l _ ^'(^k+q|| )) & (^k ^k+qy 
+ l5kq;a| 2 (l+ n ( £ k)) (l+n(£k+ q|| )) S (£k+£k-fq|| + W q) 


where the “force vertices” ffkq-n(a) re f er t° the normal 
(subscript n) and anomalous (subscript a) processes 


= («£-<+„) l k, t+qi| . q , (43) 

< q =(4t-eJ + „) (44) 


with the explicit expressions for vertices given in (B5) 
and (B6). We note that Im[M^“^(*0 + )] ^ 0 separately 
for both normal and anomalous contributions, as to be 
expected from causality. 


C. Phonon dispersion 

In the following sections we provide a discussion of 
several aspects of thermal conductivity obtained from 
the memory function approach. We are going to dis¬ 
cuss orders of magnitude estimates for the cuprates, the 
asymptotic behavior of the scattering rate, and the re¬ 
sults from the numerical evaluation of the memory func¬ 
tion. Whenever of interest, explicit reference will be 
made to parameters relevant to La 2 Cu 04 . In what fol¬ 
lows, we parametrize phonon excitations that are rele¬ 
vant for spin-phonon coupling as given by two degenerate 


modes that are polarized in-plane and have the dispersion 

Wq = ©d y) • ( 45 ) 

©£> is the Debye energy for phonons as before. Since we 
set H = ks = 1 and also assume all lengths in units of 
lattice constant unless mentioned otherwise, ©d/2 and 
the sound velocity c will be used interchangeably. 


D. Asymptotic scattering rates 


For the “normal” emission-absorption of acoustic 
phonons in Fig. [8ja) we expect four asymptotic tem¬ 
perature regimes, which are determined in part by the 
kinematics of the magnon-phonon scattering, (i) For 
T <C ©Dj thermalized magnon-phonon scattering in¬ 
volves only long-wavelength excitations, (ii) As T is 
increased, scattering from the zone-boundary, optical- 
like phonons with < T and a high density of states 
begin contribute with the Arrhenius-like, activated In¬ 
dependence. For the dispersion in (45) such phonons 
occur at the corners of the BZ at q ~ Qaf = (f, n). (iii) 
For ujq AF j(1± <C T -C J, magnons can be scattered by 
thermalized phonons from k = (0,0) to Qaf- (iv) Finally, 
one can consider a regime T J in ( |42[ ), which however 
is unphysical regarding the definition of magnons. 


1. Acoustic phonon regime 


First, we consider the “normal” scattering processes 
and r<©D. For that we may use the long-wavelength 
limit |k + qy|, |k|, |q||| •C 1 for all the in-plane wave 
vectors. As discussed in detail in Appendix [B] the in¬ 
plane modulation of J will only couple to the two out 
of three phonon branches, whose in-plane polarization 
can be always chosen as |£ q | ~ 1. We also approximate 
Wq«c|q| = C\/|q||| 2 + |qj_| 2 , with a sound velocity much 
less then the magnon velocity, i.e. c<». Expanding in 
(43) to lowest order in |k-Fq| j, |k|, and |q|| |, the memory 
function reads 


M lu ,(iQ+) = -i8 ft 


2?r (g° 


(46) 


*E 

k,q,£ 


l k ll k + qiil ( e k- e k+ q „) 




1 2 +1 q_L 1 2 

X ^(WqMek) ( 1 +^(-k+q || )) 5 (£k-£k+ q|| +Wq) 


where g^ = SX/(2y/mc) with A = dJ/dr, see Appendix | b| 
and Eq. ( |B8[ ). First, the magnon and phonon distri¬ 
bution functions n(£k) an d n(u>^) in (W6| imply that 
|k| ~ T/v and |q| ~ T/c. From this anausing energy 
conservation £k+ q|| = £k + ^q, we conclude that the in¬ 
plane phonon momentum is of magnitude |q| | ~T/v, i.e. 
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the dominant contribution of the phonon momentum is 
out-of-plane. This coincides with the conclusion reached 
after Eq. (|9j). In turn, |k + qy | ~ T/v and 
v 2 T 2 . Then, a naive power-counting, 


Sec. 


IIIB 


(J 1 _ J 1 

e k+q. 


r 


accounting also for factors of rnagnon and phonon veloci¬ 
ties, v and c, and keeping in mind that the <5(.. .)-function 
in (461 replaces one fc-integration by a factor of \/v, sug¬ 
gests that M ~ (gsp) 2 T 8 /v 6 . However, this approach 
misses one subtle detail, i.e. that the (^-integration is 
singular with 



»K) 

yiqnl 2 + |q±l 2 



which introduces an additional factor of (v /c) and finally 
leads to 



M™„(iO+) oc -iS^ ( ffs a p c ) 


p8 


(48) 


From (38) we have X/j-v ~ T 3 with a dimensionless pref¬ 
actor of order unity. Therefore, using c ~ Qd, the scat¬ 
tering rate is 


1 


oc 


Qdv 3 ’ 


(49) 


with the dimensionless prefactor (<?sp/c) 2 . This result is 
identical to our Eq. (191 and the discussion is completely 
in line with the consideration given in Sec. |IIIE| 


FIG. 9: (Color online) Difference quotient n = 

A ln[Af #H ,(iO + )]/AlnT obtained from numerical evaluation of 
A4 M „(iO + ) versus T (dots) for various Qd (lines are guides 
to the eye). Error bar is an estimated maximum “error”: 
n — 91n[M M y(iO + )]/dlnT. For reference note that for J = 
1550K, Qd = 0.1J corresponds to a maximum phonon en¬ 
ergy of approximately 270K. 


rate 


oc 


' opt 


rp 4 
w OV 2 


(51) 


2. Zone-boundary phonon regime and others 

We now consider the normal scattering processes from 
the optical-like, zone-boundary phonons with qy ss 
(■7T, 7r). First, we focus on temperatures J T wo ~ 
Qd- Calculations are simplified by shifting the zero of 
the planar component of the phonon momentum to the 
edge of the BZ, i.e. qy —> qy + (n,n), and therefore 
q —>• q + (n, n, 0). Since T< J, one can still expand in 
(43) with respect to small momenta |k + qy | and |k| near 
their respective k-points in the BZ, while for phonons one 
can set ~ ojq for all q involved in the scattering. Us¬ 
ing this expansion for vertices as discussed in Appendix 
|B| one can easily see that it leaves the structure of the 
expression for the memory function in (46) almost in¬ 


tact, with the only difference that both q and qy are 
large. Then, the power counting proceeds the same way 
as above with all the variables related to magnons gov¬ 
erned by the same smallness of the typical momentum 
k ~ T/v, while for the phonon occupation number we 
now have n(Wq) ~T/wo and the summation over q± has 
no restrictions. Altogether, this yields 


M”„(iO+)oc-iV (g£Y 


T 7 

UJ 0 V 4 


(50) 


where = 4:SX/\/2mQD as given in (B12). Similarly 


in a complete agreement with the relaxation rate ob¬ 
tained in Sec. IIIIDl Eq. |l8|). 


Without repeating similar considerations for the 
anomalous scattering in Fig. [8](b) , we simply note that 
the numerical evaluation to be presented in the next sub¬ 
section shows that for experimentally relevant tempera¬ 
tures T > 0.1 Qd its contribution is smaller by several 
orders of magnitude as compared to normal scattering. 

Finally, we note that in the unphysical regime of T^> J, 
the memory function scales as M ~ T 2 due to the three 
distribution functions and the prefactor of 1/T. Together 
with Xpv ~ T for this regime, this yields a relaxation rate 
of t -1 oc T. 


E. Numerical analysis of the memory function 

Figure [9] sheds light on our asymptotic analysis of re¬ 
laxation rates from an unbiased numerical point of view. 


The memory function has been evaluated using Eqs. (42) 


and (45) for three representative Debye energies Qd- The 


to the acoustic case, this implies the following scattering 


numerical integration has been set up to satisfy the en¬ 
ergy conserving delta-function exactly through the nu¬ 
merical solution of £k+q = £k + for each integrand 
call. This leaves a 4D integration with a non-trivial in¬ 
tegrand. Evaluation of the latter has been performed for 
the temperatures shown by the dots in Fig. [9j It shows 
a numerical approximation to (91n[M All/ (iO + )]/91nT, the 
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quantity which allows to clearly identify the temperature 
ranges where the memory matrix follows a power law, 
M fiV (iO + ) ex T n . These ranges can be seen as plateaus 
with the heights directly giving the exponent n. 

We would like to stress several points. First, as pre¬ 
dicted, a clear T 8 regime, corresponding to r” 1 ex T 5 , 
can be observed for T <C &d- However, comparing with 
the energy scales appropriate for La 2 Cu 04 , such a regime 
is unlikely to be observed in the current experimental 
studies of it, because the scattering in this temperature 
range is dominated by the grain boundaries.^ Second, 
while this may be an artifact of our scheme of modeling 
the optical phonon scattering by the phonon spectrum 
with a single phonon branch, for a robust T 7 (t" - 1 oc T 4 ) 
regime to occur, the Debye energy needs to be low enough 
compared to J. This is evident from the set of data in 
^g. m with Qjj = 0.1J where clear indications of the 
T 7 range can be observed. However, it is obvious from 
the results for <d D = 0.2J and 0£> = 0.3J, that this 
regime rapidly merges with the onset of thermalization 
of the high-energy magnons, where no well-defined expo¬ 
nent can be observed. Third, the increase of the exponent 
between the T 8 and the T 7 regimes is consistent with an 
exponential increase of M A1 „(jO + ) with temperature. This 
is exactly the signature of the Arrheni us (act ivated) be¬ 
havior mentioned as regime (ii) in Sec. IV D Finally, as 
expected for T> J we get M pv (i0 + ) <x T z . 


F. Thermal conductivity 

In this subsection we combine the memory function 
analysis of the spin-phonon coupling with the scattering 
mechanisms discussed in Sec. IfflGl that are most sig¬ 
nificant for magnon thermal transport in large-J anti- 
ferromagnets such as La 2 Cu 04 . Moreover, we will use 
parameters that may provide a reasonable description 
of the experimental thermal conductivity data, e.g. of 
Ref. [U To leading order, additional scattering mecha¬ 
nisms can be added according to the Matthiessen’s rule, 
i.e. summing their respective momentum-independent 
relaxations rates as in Sec. IIV Al 


1. Grain boundaries 


Grain boundary scattering is described in complete 
analogy with Sec. m G 1[ i.e. the results for thermal 


conductivity trivially follow from Eq. (391 with the re¬ 


laxation time r g b = lgb/v , where Z g b is the typical grain 
size and the magnon mean-free path, see (25). As was 


established in Ref. [9, this scattering dominates thermal 
conductivity of La 2 Cu 04 at T<200K and one can show 
that Eq. (39) with l g b of order of a few hundred lattice 
spacings and appropriate choice of J and lattice param¬ 
eters gives an excellent quantitative description of the 
experimental n g b(T) in this regime. For the remainder 
we choose l g b = 300 lattice constants as in Sec. Ill HI 



T [K] 

FIG. 10: (Color online) Renormalization factor of the inverse 
mean-free path due to grain boundaries from spin-phonon 
scattering for various Qo and g, J = 1550K (lines are guides 
to the eye). 


2. Grain boundaries and phonons 


We now would like to test if the magnon-phonon scat¬ 
tering can have a significant impact on the magnon mean- 
free path and thermal conductivity. For that, we calcu¬ 
late numerically the spin-phono n cu rrent relaxation rate 
(M(iO + )/x) M „ = S^/ts p from (42) as discussed above 
and neglecting anomalous scattering. Then the spin- 
phonon transport mean-free path is defined as l sp = vr sp 
and can be combined with the grain-boundary mean-free 
path Z g b to see the effect of the former. To make such 
a comparison, we also need to introduce a dimensionless 
spin-phonon coupling constant. In line with the discus¬ 
sion of the asymptotic limits in Se c. |IVD| above and in 
accord with the results of Appendix|B[ such a dimension¬ 
less coupling can be written as g=4SX/( Jy/2m&D)- The 
range of its typical values is discussed in Sec. |IIIF| and 
Appendix [C] and it is concluded that it may not exceed 
unityP 


Figure [TO] displays the ratio of the inverse magnon 
mean-free paths with and without the spin-phonon scat¬ 
tering Zgb/Ogb 1 + &(?)) f° r a ran g e °f resonable spin- 
phonon coupling constants and Debye temperatures. 
This Figure clearly demonstrates that phonons can be 
expected to be the dominant scatterers for T > 200A' 
in large-,/ antiferromagnets with a modest spin-phonon 
coupling. 


In Figure [TT] the thermal conductivity is shown for the 
scenario when only grain boundaries and phonons are 
involved in magnon scattering and for several represen¬ 
tative sets of 0£> and g. The latter are chosen to yield 
the same maximal value of k(T) at fixed / g b, with J and 
other parameters fixed to loosely match La 2 Cu 04 . Here¬ 
after all plots of the thermal conductivity display ab¬ 
solute values of n. These follow from our calculations, 
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T[K] 


FIG. 11: (Color online) Magnon thermal conductivity for 
the scattering by grain boundaries and phonons for various 
experimentally reasonable values of 0a and g, J = 1550K 
(lines are guides to the eye). 



T [K] 


FIG. 12: (Color online) Magnon thermal conductivity 

for the scattering off the effective grain boundary due to 
temperature-dependent correlation length and phonons for 
various experimentally reasonable values of Qd and g, J = 
1550K (lines are guides to the eye). 


scaling the 2D conductivity to the 3D materiaP, and us¬ 
ing IgbJks/[fiac z ) ~ 645W/(Km), for J ~ fcsl550K and 
Igb ~ 300a, with intra(inter)-plane lattice constants a(c z ) 
frorrP. We emphasize, that the magnitude of n sh own 
in Fig. |Tl] is within the typical range for La^CuO i EHiD 
The Figure shows that the spin-phonon coupling con¬ 
stants necessary to effectively suppress the conductivity 
at T > Qd are well within the acceptable values. More¬ 
over, the Figure demonstrates that a rather natural tem¬ 
perature range for the maximum conductivity to occur 
in La 2 CuC >4 and related materials is of the order of 0p. 
Comparing with Fig. [9j one can see that the decrease of 
k vs T occurs in the temperature regime where no well- 
defined power law exists for r s “ 1 (T), independently of the 
choice of the Debye energy. 

Fig. [TT| should also be compared with the upper 
curve in Fig. [6j which is obtained from the effec¬ 
tive 1/r-approximation within the Boltzmann formal¬ 
ism. Taking into account the differences between the 1/r- 
approximation in Boltzmann approach and the memory 
function calculations, different types of modeling of the 
phonon bath, and keeping in mind differences in Debye 
energies and spin-phonon coupling constants, the close 
similarity of the overall shape and other features of the 
k(T) curves from the Boltzmann and the memory func¬ 
tion approaches in Figs. [6] and m are remarkable. 


3. Grain boundaries, phonons, and finite correlation length 

Next we take into account scattering due to the finite 
spin-spin correlation length in the paramagnetic phase 
above the ordering Neel temperature, analogous to the 
discussion in Sec. [Tirol Similar to the Boltzmann ap¬ 
proach, a formally proper treatment of its effects is be¬ 


yond the present study. A qualitative description, how¬ 
ever, can be obtained readily. First, a finite correlation 
length implies a ma ss gap in the magnon dispersion, dis- 
see Eq. (28) Second, since the no¬ 


cussed in Sec. Ill G 


tion of magnetic order is meaningful only within patches 
of linear dimension ~ £(T), Eq. (26), the system con¬ 
sists of effective “grains” with a temperature-dependent 
size l g b{T) « min(^ g b, £(T)), the sentiment expressed ear¬ 
lier in Sec. |III| Once again, Matthiessen’s rule is used to 
yield a smooth crossover between the two grain-boundary 
regimes as a function of temperature and to include the 
spin-phonon scattering with the inverse mean-free path 
given by l~\T) = l~£ + ^\T) + l~ p \ 

First, recalculating n within the relaxation time ap¬ 
proximation, Eq. (39), with a constant scattering time r 
and the mass gap due to correlation-length in the magnon 
dispersion, Eq. (281, affects k only weakly and only at 
higher temperatures, in a broad agreement with the con¬ 
clusions reached in Sec. |III H| For brevity, we will not 
study the impact of the mass gap on the memory func¬ 
tion. In turn, the main effect of the correlation length 
on the thermal conductivity is from the limiting of the 
mean-free path by £(T). For the parameters relevant to 
the cuprates, the correlation length gets short enough to 
likely dominate dominate the thermal conductivity for 
T > 250K. 

Our Figures [l2| and [Ts] display the combined effect of 
spin-phonon and effective boundary scattering, discard¬ 
ing the magnon mass gap. Fig. [12] shows that at high 
temperatures, in contrast to the rather slow decrease of 
thermal conductivity in Fig. 11 k(T) is dominated by 
the exponential decrease of the effective grain size £(T). 
This figure also shows that small differences in the re¬ 
sults Fig. [TT] due to differences in coupling constants and 
Debye energies are completely suppressed by inclusion 
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T [K] 

FIG. 13: (Color online) Contrasting impact of different scat¬ 
tering mechanisms on the thermal conductivity. Phonons and 
crystalline grain boundary (upper curve), correlation length 
and crystalline grain boundary (middle curve), phonons, cor¬ 
relation length, and crystalline grain boundary (lower curve). 


of £(T). Fig. 13 contrasts the various scattering mech¬ 
anisms against each other. This figure clearly demon¬ 
strates that while the spin-spin correlation length seems 
to provide the major limit on k(T ) at higher tempera¬ 
tures, magnon-phonon scattering still contributes signifi¬ 
cantly to the suppression of the magnon heat current up 
to high temperatures. 

Once more we emphasize the close similarity of the 
combined results from the memory function method in 
Fig. [13] with those from the Boltzmann theory in Fig. [7] 
Given that these results have been obtained indepen¬ 
dently from two distinct theoretical approaches, this 
agreement strongly corroborates our main conclusions. 


V. SUMMARY 

To summarize, we have considered relaxation mecha¬ 
nisms of 2D magnons in large- J antiferromagnets such as 
La 2 Cu 04 with the goal of providing a basis for quanti¬ 
tative calculations of thermal transport by spin degrees 
of freedom in this class of materials. We conclude that 
the magnon thermal conductivity in these and related 
systems is limited by three key scattering mechanisms: 
grain boundaries, 2D correlation length, and magnon- 
phonon scattering, the latter being effective at interme¬ 
diate and high temperatures, T > Op. Impurity-like 
lattice-induced disorder and magnon Umklapp scatter¬ 
ing have been found to be less important. 

The bulk of our consideration has been devoted to the 
scattering of 2D magnons on 3D phonons, with acoustic 
and optical phonon branches examined in detail. Within 
the Boltzmann approach and 1/r approximation we have 
advocated the use of a simplified “effective” phonon DoS 
approach, which allows for straightforward yet fairly re¬ 


alistic calculations yielding an effective expression for the 
magnon relaxation rate on optical and acoustic phonons 
that contains minimal number of parameters. We have 
also employed the memory function approach, in which 
we have retained the full microscopic expressions for the 
magnon and phonon spectra, with the only approxima¬ 
tion being the coupling of spins to only a single dispersive 
phonon branch. 

Within both approaches, we have closely analyzed the 
power-law asymptotic regimes of the magnon relaxation 
rates involving acoustic and optical, or zone-boundary, 
phonons. Not only have such investigations proven very 
instructive in the analysis of which part of the magnon 
population is affected most strongly by phonons, but 
also have demonstrated a remarkable accord between the 
two very different theoretical approaches to the transport 
problem. 

In both Boltzmann and memory functions approaches, 
we have demonstrated that having a modest spin-phonon 
coupling, well within a reasonable range of parameters, 
phonons can have a significant impact on the magnon 
mean-free path and thermal conductivity. Taking into 
account other scattering mechanisms within the ther¬ 
mal conductivity calculations, we have concluded that 
magnon thermal conductivity in large-J antiferromag¬ 
nets should be largely controlled by the “boundary-like” 
scatterings, coming from either the real grain boundaries 
or the 2D correlation length, with a substantial correc¬ 
tion from the magnon-phonon scattering, affecting n(T) 
at intermediate and high temperatures, T>Qjj. We have 
also demonstrated that the results of both approaches on 
thermal conductivity are remarkably similar. 
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Appendix A: Spin model 

For the spin system we consider the two-dimensional, 
nearest-neighbor, spin-1/2 Heisenberg antiferromagnet 
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on a square-lattice 


^ s = J^S i -S i , (Al) 

(ij) 


where summation is over the nearest-neighbor bonds and 
J is the exchange coupling constant. We treat this model 
using standard linear spin-wave theory. For the bipartite 
lattice we transform from the laboratory to the rotated 
frame for spins and bosonize spin operators according to 

S-° = e iC * Ti S* , S*° = e lQr 'S*, S?° = Sf, (A2) 
S? = S- b\bi, S+ = (S~y » biV2S, 

where Q = (n,ir) is the antiferromagnetic ordering vec¬ 
tor. Then, within the 1/S approximation, one needs to 
retain only quadratic terms in the Hamiltonian © 

Si ■ S l+5 => S[b\bi + bl +s b l+5 - b\ +5 b\ - b i+5 b t }. (A3) 


After the Fourier transform, the subsequent treatment 
of (Al) involves Bogolyubov transformation of magnons 
given by 


where A = adJ/dr, with a being the lattice constant, <5 
runs over the nearest neighbors, S are the correspond¬ 
ing unit vectors of the lattice, AU i]( 5 = (U i+ j—Uj), and 
U, is the displacement vector of the ith ion. While the 
Hamiltonian in ( |B1[ ) is rather general, one can see it as 
describing a single 2D plane of spins on a square lattice, 
representative of a Cu0 2 plane of La 2 Cu 04 . Below, we 
will extend this picture to a tetragonal array of magneti¬ 
cally decoupled planes to take into account the 3D nature 
of phonons. Note that while t he d isplacement of ions are 
described by 3D vectors, in (Bll only their projections 


on the in-plane bonds (5) matter. 


After bosonizing spin operators according to (A3), us¬ 


ing the symmetry of the square lattice, Fourier transform, 
and some algebra, we obtain 


H s . ph = HAS E £ (* • IJ. 

q.k.k' J 

k-k’ 

A1 A 


COS 


1 


-COS 

2 


2 

k + k' 


I Sill 


k k' 


s ) &L&k 


(B2) 


(^-k + fc-k'&k) 


bk=u k /3 k +v k ^_ k , (A4) 


with 


v-k 


1 + Dc 
2i/ k 


u k = sign (q k ) 




1 ~ »k 

2vk 


(AS) 


where jz k = yH — is related to the magnon energy via 
e k = 4JSV k and 7 k = (cos fc^+cos k y ) /2. Finally, the 
spin-wave Hamiltonian is given by 

=J2e k f3l(3 w , (A6) 

k 

At long wavelength « vk with the spin-wave velocity 
v = V2J. 


Appendix B: Microscopic derivation of spin-phonon 
Hamiltonians: vertices, polarizations, etc. 

a. Spin-phonon Hamiltonian: Bravias lattice 


where q = (k — k',(jj_) and the summation over the di¬ 
rections of the bond now takes the values of only two 
unit vectors, S = x and y. I n derivation of (B2) we have 


performed summation over the sites of the tetragonal lat¬ 
tice of magnetically decoupled planes. It is natural that 
the in-plane component of the phonon momentum (in the 
index of U q ) is tied to the magnon momenta via momen¬ 
tum conservation, q | = k — k', while the component of 
the phonon momentum perpendicular to the plane, qu, is 
not conserved and is an independent variable in the sum 
in (B2|. This is because the 2D magnons have zero cou¬ 


pling between the planes (infinite mass in that direction) 
and thus provide no restriction on q^. 

The subsequent treatment of ( |B2| involves Bogolyubov 
transformation of magnons as in (A41 and rewriting the 


displacement operators in terms of phonons. For the 
Bravias lattice, the q’s Fourier component of the lattice 
displacement is given by 


u q = E 




2niWq 


i (B3) 


For an isotropic, nearest-neighbor interaction of spins 
and considering the lattice of only magnetic ions, the 
most natural source of the spin-phonon coupling is from 
the modifications of the superexchange due to local 
stretching or compression of the bond length^ 


Hs-ph = E SJ ( r i ~ r ») S * ‘ S i 
(ij) 

» \ E ( s ■ au *a) s, • S i+4 , 

i : 5 


where £ = 1,2,3 numerates one longitudinal and two 
transverse phonon polarizations, are the polarization 
unit vectors, m is the mass of the unit cell, and w q and 
u q £ are the energies and operators of the corresponding 
phonon branches, respectively. 

Applying (|B3[) and (|A5|) to (B2) we obtain 


n s - ph = E Ef^kcq/^k 

k, k ',qx £ 


(B4) 


+ 2 Kk' q (&0lk + H.C.) } (« q t + «-«*) 


(Bl) 
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where “normal” and “anomalous” magnon-phonon ver- 


(B5) 


tices are 




v e 

^k,k',q — 

4SA 

cz=x,y 

q~ 

^2 moo ^ 

2 

X < 

f q~ 
l C0S T 

(uu -(- vv') — 

COS 

r rOd,^ 

^k,k',q — 

4SX 

E #“ si * 

ce=x,y 

q~ 

yjzmuo^ 

2 

X < 

f q~ 
l C0S T 

(uv' + vu') — 

COS 


(B6) 


where we have absorbed i in the definition of phonon 
operators, and are the x and y projections of the 
corresponding vectors, and we have introduced shorthand 
notations q±=k±k', =u k o), and v^=v k i>). 


b. Asymptotics and polarizations 

Of interest, of course, is the asymptotic behavior of 
these vertices in the vicinity of k, k' —► T, Q as the mo¬ 
menta of magnons are confined to these regions at not 
too high temp erat ures. 

Expanding (B5) for magnon momenta k. k' —> T (or 
Q), which also corresponds to q|| —>• 0, gives a = x,y 
components of the vertex 

SX^ a 


v. 


k,k',q 


y // 4 


1 + 


2 k a k' a 

w 


(B7) 


7-od,£,ct 


ft,a. 


and expansion of (B6) yields 1j] 

In general, three polarization vectors £q in (|B3|) 


„ = —V 

k.k',q k.k'.q’ 


should be obtained as solutions of the dynamical ma¬ 
trix equation for the lattice vibrations of the tetragonal 
Bravias lattice for each q-vector. However, in the long- 
wavelength limit q—>-0, these solutions can be classified 
as in the continuum as one longitudinal and two trans¬ 
verse modes. For the former, £q ?) ft q/|q| is along the 
momentum of the phonon, while the latter can be chosen 
freely as an orthogonal pair in the plane perpendicular to 
q. For instance, a convenient choice is £q * ft y x q/|y x q|, 
which ensures that it is orthogonal to the y — q plane. 
Since in (B5) (B7) we are interested only in x and y 


projections of the polarization vectors, this choice leaves 
nonzero only component. 

However, the situation is even simpler because, accord¬ 
ing to the discussion in Sec. |IIIB[ in a typical scattering 
process component of the phonon momentum perpen¬ 
dicular to the plane is much larger than the in-plane 
component, |qj_| |q|||- This means that the in-plane 

projections of the longitudinal phonon polarization vec¬ 
tor are negligible, ~q x ^/q± <C 1. In turn, this 

means that the polarizations of the two transverse modes 
lie almost entirely in the x — y plane and the correspond¬ 
ing vectors can simply be chosen along the x and y axes 


with ^q^’ 


, A2),v, 

' Sq 


: 1 and ^ 1),y * 


, f(2)v 

' Cq 


0. Therefore, 

it is the transverse phonons with the momentum largely 
orthogonal to the antiferromagnetic planes that couple 
most strongly to the spin excitations. As a result, sum¬ 
mation over l in (B4) should concern only them. We note 


that the same arguments have been used previously for 
the spin-phonon coupling in ID spin chains, see Ref. 1171 


Thus, using uo, 


(i) 


oo, 


( 2 ) 


c|q| for the transverse 


branches, the magnon-phonon coupling with the acoustic 
phonons is reduced to just two non-zero terms 


Tyl 1 ).^ „ac ^11 

W.k'.q ~ Osp 


V < ' 2 ' > ’ V ft O aC 

^k.k'.q "sp 


y 


, _ . 9 h x 1c ,x \ 

v»(i + w)’ ( B8 > 


% 


/fq| 


VWm (i + 


2 k y k ,y \ 


where g “ = SX/2^/mc and the results are the same for 
the anomalous coupling V kk ! . 


Since the vertices in (B 81 couple to different but de¬ 


generate phonon branches, they will contribute indepen¬ 
dently to the scattering probability, so that |Vtot | 2 = 
lyl 1 ),®! 2 _|_ |pr( 2 ),j/| 2 Then p; j s natural to introduce an 
effective coupling to just one phonon branch that would 
lead to the same scattering rate 


Vi 


eff 

k.k', 


9s 


x/fqf 




(B9) 


in which we simply ignored the angular dependence in the 
brackets in ( |B8| [second terms]. This is equivalent only to 
a quantitative (order of 1) change in the e ffect ive coupling 
constant. Needless to say, the result in (B9) is identical 
to the coupling proposed in Q and used throughout the 
paper. 

With that the effective Hamiltonian for magnon- 
phonon coupling becomes 


nf_ ph = E ^,q{/ ? ^ k + 2 (^^ k + H -c-)} 

X + «-q t) > (BIO) 


k,k',(j ± 


c. Zone-boundary phonon 

A yet another asymptotic consideration is relevant for 
the magnon-phonon coupling on a Bravias lattice. It con¬ 
cerns scattering of a magnon between two branches, from 
the vicinity of T point to the gapless branch at the anti¬ 
ferromagnetic ordering vector Q or vice versa. Therefore, 
of interest are the limits k—and k'->Q for vertices in 
(B5) and (B6). Importantly, the phonon involved in such 


a process has a large momentum with the in-plane com¬ 
ponent q|| ft— Q, which corresponds to a zone-boundary 
excitation with the energy of order ©d and is similar to 
the optical modes considered next. 
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Expansion of the “magnon part” [content of the curly 
brackets] in (B5) and (B 6 ) gives 


W.k'.q 


T7-od,£ 

*K,k',q 


4SA { k '<*\k\ - k a \k'\ 

a q l 2^/fkljk 7 ] 

4^A ^^ a { k a \k'\ + k' a \k\ 

q i avikiM 


(BH) 


where we have shifted k' —> Q + k' and also used that 
sin^y*/2) ~ —1 for qy sa— Q. Focusing now on the “nor¬ 
mal” vertex and following the same choice of polarization 
vectors as above leads to two non-zero vertices 


T/(!)»* r- 

K k,k',q ^ 

» q zh 
y sp 

t /(2),2/ 

K k,k',q ^ 

» q zh 
ysp 


k ,x |k| - k x |k' 


2^/\kW 


(B12) 


with the coupling constant to the zone-boundary phonon 
g z p = 4S'A/V2 to©d- Combining the two couplings to 
different phonon branches into one effective vertex gives 


4k.k',q ~ 9s p V\km sin (^/2). 


(B13) 


where ip is the angle between k and k' and anomalous 
vertex is the same with sin— > cos. The only difference of 
the effective coupling in (B13l from the coupling to the 
optical phonon proposed in (|7j is the additional angular 
dependence. It is easy to see that when calculating the 
scattering probabilities with |V | 2 the latter will be aver¬ 
aged to an additional factor 1/2 and thus corresponds to 
a simple change of the effective coupling constant. 


d. Optical phonons 


To introduce optical phonons into the spin-lattice cou¬ 
pling model one needs to depart from the Bravias pic¬ 
ture of (Bl). For that, keeping in mind cuprates, mu¬ 


tual displacement of copper ions AU i>( 5 = (Ui + ^ — U,) 
should now be recognized as involving the complete set 
of phonon normal modes. That is, the summation over 


i in (B31 should now be treated as involving not only 


three acoustic branches of the Bravias lattice, but the 
entire group of optical modes as well. In that sense, the 
general form of the spin-coupling Hamiltonian in (B2) 


and the entire consideration of it remains valid with the 
magnon-phonon coupling constants and optical phonon 
energies used accordingly. 

It is then clear that the small-qy scatterings involv¬ 
ing optical branches should be less important than the 
ones involving acoustic modes, while the large-qy do not 
carry additional smallness of |qn | ~ T/v and should be im¬ 
portant. Therefore, a consideration of the leading effect 
of the optical-phonon coupling is very similar to that of 
the zone-boundary phonon above. Thus, without going 
into consideration of the details of the crystal structure 


of specific materials and neglecting nonessential angular 
dependence similarly to the cases considered above yields 

^k,k', q » 9T VWM, (B14) 

with g° -P* S3 4\S/y/2muj 0 where w 0 is the energy of the 
optical branch. Thus, the form of the magnon-phonon 
interaction proposed in ([7]) is verified. While, obviously, 
the coupling strength must depend on the type of the 
optical model involved, one can suggest the so-called 
“stretching mode” at high energies as a strong candidate 
for a significant spin-phonon coupling. 


Appendix C: Magnitude of spin phonon coupling 


All three effective magnon-phonon coupling constants, 
to acoustic, to zon e-bo u ndary , and t o opti cal modes, in¬ 
troduced in Eqs. (B 8 ), (B12|, and (B14), respectively, 
have very similar structure: g sp oc A/ y/2 uiQd with a co¬ 
efficient of order of unity. Here A = adJ/dr ss yj is the 
response of the superexchange constant to the atomic 
displacement. It has been arguecP^ that because the su¬ 
perexchange is very sensitive to the interatomic distance, 
the typical values of 7 are ~ 10 — 20. However, this large¬ 
ness is offset by the smallness of a characteristic scale as¬ 
sociated with phonons,^ l/y/m&D^ 1/100. This is, ac¬ 
tually, the same parameter that characterizes the small¬ 
ness of the typical magnitude of the zero-point atomic 
displacement relative to the interatomic distance^ or the 
smallness of the typical velocity of an atom in a lattice 
relative to the sound velocity. Thus, the physical range 
of the magnon-phonon constants is g sp /J~ 0 . 1 . 

To make a closer estimate of the spin-phonon coupling 
constants as related to the cuprates, we use m = m c u 
and a typical phonon Debye energy Q D = 400K. The 
largest uncert ainty is in the value of dJ/dr = 7 J /a with 
0(3... U'pM!. Therefore 


4 SXh 

y/2ma 2 kBQD 


= gJ ~ 0(0.05... 0.3) J, 


(Cl) 


where we have used a ~ 3.8A, as in I^CuCH. This 
defines the dimensionless spin-phonon coupling constant 
g and demonstrates that spin-phonon coupling in the 
cuprates, while significant, is still within the bounds to 
justify the use of perturbation theory. 


Appendix D: Verification of the “effective phonon 
DoS” approach 


Here we complement the discussion of Secs. and 
|III E| by providing more verifications of the effective DoS 
approach advocated in this work. 
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FIG. 14: (Color online) Magnon relaxation rate on the optical 
phonons, Eq. (161, vs £k for a representative T = 400K and for 


effective phonon DoS Models II and III in Eqs. (13 ) and (141. 


The results are normalized to the high-temperature asymp¬ 
totic behavior of the phonon-absorption term t^ 2 \ Eq. JlTp, 
[T 2 s£/v 2 wo]. Parameters are indicated in the graph. The 
vertical axis is in units of (g°p t /' y ) 2 - Dashed line is the result 
of the direct numerical integration in Eq. § for the disper¬ 
sive optical phonon in Fig. [2] without the approximation of 
Eq. (101. Inset shows individual contributions of the three 
diagrams in Fig. [I] for both the effective DoS approach (161 
with the Model II and the direct numerical integration. 

e. Optical phonons 


For the “Model III” in Eq. (14) (flat phonon mode) 


of the effective phonon DoS the integral in Eq. (16) is 


trivially removed and the relaxation rate is given by a 
compact analytical expression 


1 

7k 

x 


° P t 

i/sp 


{©(, 


£k ^ i 


( 5 ) < D1 > 

w 0 ) (£k - w 0 ) 2 (n{uo) + n(s k - w 0 ) + l) 

+ (£k + wo ) 2 ^n(wo) — n{e k + wo)^ 

+ 0(w o - £k) (w 0 - £k) 2 (n{u 0 - £k) - ra(w 0 )) }. 


An identical expression can be obtained directly from 
Eq. ([5 ) for the optical phonon energy w q = Wo and lin- 
earizec magnon energy and magnon-phonon vertex in 0 
and 0. Thus, in this case, the effective phonon DoS ap¬ 
proach is exact. 

Fig. 0 shows the expected activated and high- 
temperature asymptotic behavior of the scattering rate 
vs T from fl8| ). Fig. [14] demonstrates the validity of an¬ 
other aspect of the asymptotic consideration in Eq. ( fl7] ) : 
the dependence of 1/Vk on £k at T>wo- The results in 
Fig. 14 are normalized to the high-temperature asymp- 


( 2") 

totic behavior of the phonon-absorption term I/t^ in 


Eq. (|17|), T 2 £^/v 2 u} 0 . Clearly, for £ k <T such a behavior 


FIG. 15: (Color online) T-dependence of the magnon relax¬ 
ation rate on the acoustic phonons, Eq. (211, for £k = T and 
using effective phonon DoS model in Eq. G0-, The results are 
normalized to the asymptotic behavior, Eq. (191, [T 5 /d 3 wo]. 


Parameters are as discussed in text, Debye energy @£> =400K. 
The vertical axis is in units of (<?J p/v) 2 . Individual contribu¬ 
tions of the first two terms in (21 1 [diagrams in Fig. 0a) and 
(b)] are indicated. Dashed lines are the result of the direct 
numerical 3D integration in Eq. ( |5|) fo r the acoustic phonon 
without the approximation of Eq. (10[ ). 

the phonon-emission term, 1 /t^\ carries higher power 
of £k. also in agreement with Eq. In addition, the 

results of a direct 3D numerical integration in Eq. 0 for 
the dispersive optical phonon in Fig. [2j w q w wo + aq 2 , 
without the approximation of Eq. ( |10[ ) are shown by the 
dashed line. Inset shows individual contributions of three 
diagrams in Fig. 0 for both the effective DoS approach 
(161 with the Model II and the direct numerical integra¬ 
tion. One can see a very close agreement of the effec¬ 
tive DoS method with the direct numerical integration 
in 0 , which is achieved at a fraction of numerical cost 
as the former approach requires only a ID integration in 
Eq. @. 


/. Acoustic phonons 

With the help of our Figs. [0] and [0] we provide a 
demonstration of the accuracy and numerical efficiency 
of the effective phonon DoS approach. In them the T- 
and £k-dependencies of the relaxation rate by acoustic 
phonons from (21) are compared with an explicit numer¬ 
ical 3D integration in 0 with the magnon-phonon ver¬ 
tex from ® and without the use of approximation (10). 


These figures also offer an additional confirmation of the 
asymptotic trends of Eqs. (19) and (22). 

In Fig. [15] the T-dependence is shown for the relax¬ 


ation rate Eq. (21) for £k = T, i.e., on “thermal shell”. 


The results are normalized to T 5 /v 3 Qd to make the 
asymptotic behavior of (19) apparent. The vertical axis 
is in units of {g^/v) 2 . We also show individual contri- 


is confirmed (finite intersect of the vertical axis), while butions of the first two terms in (21), while the third 
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FIG. 16: (Color online) Magnon relaxation rate, Eq. ( |21| ), 
vs £k for a representative T = 400K and for effective phonon 
DoS model in Eq. ( |15| ) . The results are normalized to the 
asymptotic behavior, Eq. (22), [T£^/v 3 Qd]- Inset: Same with 
the normalization by T 2 £^/v 3 Od- 


is negligible as discussed above. The results of a direct 
3D numerical integration in Eq. (I5| for acoustic phonon 
(within the Debye approximation) without the approxi¬ 


mation of Eq. (101 are shown by dashed lines. One can 
see a very close quantitative agreement of the effective 
DoS method with the direct numerical integration. At 
very low T, the direct numerical procedure becomes un¬ 
reliable due to very small k'- and q-space of integration 
relevant for the scattering. Again, the high-accuracy re¬ 
sults of the effective DoS approach are achieved at a frac¬ 
tion of the numerical cost of the direct integration. 


Fig. 16 shows the dependence of 1/Vk on £k at a rep¬ 
resentative T = 400K. The results in Fi g. [I6l are normal¬ 
ized to the asymptotic behavior in Eq. (|22[), Te^/v^Qd- 
Clearly, for 1 /t^ such an asymptote is essentially pre¬ 
cise for all the energies, while for the phonon-absorption 
term substantial deviations occur at lower energies, the 
feature also emphasized in the inset of Fig. [16] This can 
be analyzed by a more careful asymptotic treatment of 
Eq. (21) in the £kregime, which show that the small¬ 
ness of the subleading terms by c/v gets compensated by 
the largeness of T/e k for small enough £k, and the ul¬ 
timate asymptotic behavior of this te rm in th e £k -C T 


regime is 1 /t( 1 2 ) ~T 3 £^./u 4 5 6 * 8 9 10 11 12 . As in Fig. 
a very close quantitative agreement of 


15 


Fig. 


16 


shows 


the effective DoS 
method with the direct numerical integration. 
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